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ABSTRACT 

The observed dynamics of gas and stars on galactic and larger scales cannot be accounted for by 
self-gravity, indicating that there are large quantities of unseen matter, or that gravity is non-Newtonian 
in these regimes. Milgrom's MOdified Newtonian Dynamics (MOND) postulates that Newton's laws are 
modified at very low acceleration, and can account for the rotation curves of galaxies and some other 
astrophysical observations, without dark matter. Here we apply MOND to two independent physical 
systems: Lya absorbers and galaxy clusters. While physically distinct, both are simple hydrodynamical 
systems with characteristic accelerations in the MOND regime. We find that, because MOND violates 
the strong equivalance principle, the properties of Lya absorbers depend strongly upon the (unknown) 
background acceleration field in which they are embedded. If this field is small compared to their internal 
accelerations, then the absorbers are more dense and about ten times smaller than in Newtonian gravity 
with dark matter, in conflict with sizes inferred from quasar pair studies. If, however, the background 
field is rather large, then the absorbers take on properties similar to those predicted in the CDM picture. 
In clusters MOND appears to explain the observed (baryonic) mass-temperature relation. However, given 
observed gas density and enclosed mass profiles and the assumption of hydrostatic equilibrium, MOND 
predicts radial temperature profiles which disagree badly with observations. We show this explicitly for 
the Virgo, Abell 2199 and Coma clusters, but the results are general, and seem very difficult to avoid. 
If this discrepancy is to be resolved by positing additional (presumably baryonic) dark matter, then this 
dark matter must have ^1 — 3 times the cluster gas mass within 1 Mpc, and about ten times the gas 
mass with 200 kpc. This result strongly disfavors MOND as an alternative to dark matter. 

Subject headings: cosmology: theory - gravitation - dark matter - galaxies: clusters: general - 
intergalactic medium hydrodynamics 



1. INTRODUCTION 

The currently most widely-accepted 'standard model' of 
cosmology holds that the vast majority of the mass den- 
sity of the universe is hidden in dark forms. The rota- 
tion curves of galaxies and the dynamics of galaxy clusters 
cannot be accounted for by the gravitation of visible stars 
and gas, while constraints from primordial nucleosynthesis 
studies imply that the additional 'dark matter' postulated 
to remedy this discrepancy must be non-baryonic. On a 
cosmological level, collisionless (or very weakly collisional) 
dark matter is required if primordial density perturbations 
of amplitude Ap/ p ~ 10 -5 are to grow quickly enough to 
form galaxies by the present epoch. Finally, recent deter- 
minations of the high redshift type la supernova Hubble 
diagram, in tandem with microwave background data im- 
plying a flat cosmic geometry, and a large collection of data 
indicating that clustering matter only contributes ~ 30% 
of the critical density, imply that the universe also contains 
'dark energy' of a yet more exotic form which causes ac- 
celeration in the cosmic expansion (see, e.g., Peebles 1999 
and Turner 1999 for recent reviews). 

Discomfort with this repeated postulation of invisible 
matter with increasingly unusual properties has led some, 
quite reasonably, to ask whether the observed phenom- 
ena could be accounted for not by the presence of unseen 
matter, but by a departure from Newtonian/Einsteinian 
dynamics in the regime where dark matter is hypothe- 



sized to be important. Perhaps the most successful of such 
proposals is the MOdified Newtonian Dynamics (MOND) 
proposed by Milgrom (1983a, b,c), in which Newtonian dy- 
namics breaks down below an acceleration threshold of 
ao ~ 10 -8 cms -2 . Successes of this theory are that it 
accounts for the rotation curves of galaxies of various lu- 
minosities (Sanders & Verheijen 1998) and surface bright- 
nesses (de Blok k McGaugh 1998; McGaugh & de Blok 
1998), accounts for the Tully-Fisher and Faber- Jackson re- 
lations (van den Bosch & Dalcanton 2000; Sanders 1996, 
2000), and - arguably - roughly accounts for the amount 
of dark matter inferred in clusters (Sanders 1999). 

Despite some attempts, MOND has not been general- 
ized into a satisfactory relativistic theory which can yield 
unambiguous cosmological predictions. Scott et al. (2001) 
explore the difficulties in applying MOND to cosmology 
and review a number of claimed conceptual and empirical 
difficulties with MOND. But in light of MOND's general 
success when applied to galaxies, and the current lack of 
any decisive empirical argument against it, it is worth in- 
vestigating whether MOND works in detail in systems for 
which it was not designed yet makes relatively unambigu- 
ous predictions, and for which good observations are avail- 
able. We propose and perform two such tests. First, we 
deduce basic properties of the Lya absorbers, making use 
of the technique of Schaye (2001). These can be compared 
directly to observations concerning the sizes and number 
densities of the absorbers. Second, we derive relations be- 
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tween the density, enclosed mass, and temperature profiles 
of galaxy clusters which can be directly compared to avail- 
able X-ray data. Both Lya absorbers and clusters con- 
stitute relatively simple physical systems, are well within 
the MOND regime, and (in the absence of dark matter) 
are dominated by gas which can be accurately observed, 
making them ideal testing grounds for MOND. 

2. MOND AND THE LYa FOREST 

As argued by Schaye (2001; hereafter SOI), basic prop- 
erties of the gas responsible for Lya absorption in quasar 
spectra can be deduced using simple physical arguments. 
Along any sightlinc passing through a gas 'cloud', the size 
of the region dominating the absorption will typically be 
of order the local Jeans length, regardless of the overall 
shape of the cloud and regardless of whether the cloud as 
a whole is in dynamical equilibrium. Using this reasoning, 
SOI calculates properties of the Lya absorbers which are in 
good agreement with numerical simulations and available 
observations. Because the absorbers are dynamically sim- 
ple and have very low characteristic accelerations, MOND 
makes strong predictions about their properties, which we 
will now derive in parallel to the treatment by SOI. 

MOND can be formulated in a number of ways, either as 
a modification of inertia (e.g., Milgrom 1999) or of gravity 
(e.g., Bekenstein & Milgrom 1984). In its most general 
incarnation, the gravitational acceleration in an 'isolated' 
(meaning not embedded in a gravitational field with larger 
characteristic acceleration) system is given by 



a = ^a N a 



(1) 



when <C an, where a at is the acceleration calculated 
using Newtonian gravity, and an ~ 1.2 x 10~ 8 cms -2 (Mc- 
Gaugh & de Blok 1998) is the MOND acceleration param- 
eter. (Systems which are embedded in a large external 
field are discussed below). This acceleration law yields a 
dynamical time 
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where L and p are the characteristic size and density of the 
system, respectively, ur is the hydrogen number density 
in a medium of hydrogen mass fraction (1 — Y), and m p 
is the proton mass. The system's sound-crossing time is 
unchanged by MOND and is 
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where T is the temperature, 7 = 5/3 is the assumed adia- 
batic index and fi w 0.59 is the mean molecular weight 
per particle for a fully ionized primordial plasma with 
Y = 0.24. Setting these timescales equal to each other 
yields the Jean length 
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This can be converted to a column density Njjj — LjUh x 
(niii/nii) using the ionization correction from SOI, valid 

1 We assume that whatever cosmology MOND engenders will be a 
cosmology with S7 m = 1 — = 0-3 (as indicated by observations). 



for a highly ionized, optically thin plasma: 
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where (3 w 4 x 10- 13 7y - 76 cm 3 s" 1 and Y = T 12 x 10~ 12 s 
are the recombination and ionization rates with T4 = 
T/10 4 K; Ti2 w 1 is measured at redshift z ~ 3 (see Scott 
et al. 2000 and references therein). The resulting Jeans 
column density can be expressed as a function of overden- 
sity S = nn/nu — 1, as 
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where T4 = T/10 K, p\ is the current critical density and 
fib is the baryonic density parameter. Thus isolated Lya 
absorbers of density equal to the cosmic mean have about 
ten times lower column density than in the CDM picture, 
with different dependences on T, z, etc. (c.f. S01, Eq. 
10). One can also express the Jeans length in terms of the 
observed HI column density: 
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about a factor of ten smaller than an absorber of the same 
column density in the CDM picture (S01, Eq. 12), and 
with different scalings for all parameters. Note that sig- 
nificant external pressure would only decrease this size. A 
few self-consistency checks are in order. First, the Newto- 
nian acceleration is 

a N w 4 x 10- 13 cms- 2 (A H/ /10 14 c m - 2 ) 2/5 T 4 °- 70 r 2 / 5 , 

(8) 

verifying that the system is in the MOND regime. The in- 
ternal acceleration is then w 00/ 170. Second, the system's 
dynamical time is 



t dyn = 2 x 10 16 s (AWl0 14 cm- 2 ) ^ T^T^'\ 



(9) 



whereas the Hubble time 1 is w (1.1 — 2.7)/i 65 1 x 10 17 s for 
3 ^ z ^ 1, so the absorbers are self-consistently in lo- 
cal hydrostatic equilibrium (assuming they can reach this 
equilibrium; see below). 

Fricdmann model with a scale factor that evolves roughly as in a standard 



The method of SOI can also be used to 'invert' the ob- 
served column density distribution into an estimate of the 
total mass density in Lya absorbers, using 

8nGH(z) ( "fkTrripT \ 2/5 
gas " 3tf 2 c(l + zf \^1-Y)(l-Y/2)P) 

^-^j^f^l (1 0) 

where d 2 n(NHi, z) / dNnidz is the differential number of 
observed lines of column density Nhi at redshift z. 

The preceding calculation applies to an isolated ab- 
sorber, but because the strong equivalence principle is vi- 
olated in MOND, the internal dynamics of systems can 
change if they are embedded in an external acceleration 
field go, even if it is homogeneous (i.e., go = goz with co- 
ordinates such that z is the external field direction). Mil- 
grom (1986) shows that in the Lagrangian formulation of 
MOND (Bekenstein & Milgrom 1984), if g is small com- 
pared to ao but large compared to the internal accelera- 
tions of a system, then gravity in the subsystem is Newto- 
nian in the coordinate system {x 1 ', y', z'} = {x, y, z/2} with 
an effective Newton's constant G' = (ao/go)G. In this 
case, properties of the absorbers in MOND can be com- 
puted directly from the corresponding Newtonian/CDM 
expressions from SOI, with those substitutions. 2 For ex- 
ample, the size of the absorbing region is then 

L 7 ~40kpc^ i0lW2 j T 12 T 4 

in the direction perpendicular to the external field, and 
somewhat (about two times) longer in the parallel direc- 
tion. The density parameter in Lya absorption systems is 
given by the CDM calculation of SOI, adjusted by a factor 

of (ao/.9o)~ 1//3 /g 1 ^ 3 , where f g « 0.16 is the gas-to-matter 
mass ratio in the CDM calculation. 

The above analysis gives two quite testable predictions 
concerning the Lya forest in MOND. First, the total den- 
sity of gas can be computed using Eq. pli (for the isolated 
case) or Eq. 16 of SOI (for the extcrnaTficld case). Since 
Lya systems are very deep in the MOND regime, there 
is no a priori reason to expect that this density will be 
at all reasonable (i.e., as compared to the nucleosynthesis 
value of fif,), and this can be assessed. Second, the char- 
acteristic sizes of absorbers, given by Eqs. ^ and O, can 
be compared to observations of lensed quasars and quasar 
pairs that constrain the transverse sizes of absorbers. In 
making both comparisons, the external field case requires 
a value for the mean acceleration field of Lya absorbers at 
z — 3. This is currently not calculable, but may be crudely 
estimated. At z — 0, typical mean observed accelerations 
can be obtained by dividing typical bulk flow velocities 
of ~ 600 km s^ 1 (see, e.g., Dekel et al. 1999; Dale et al. 
1999) by a Hubble time, yielding go/ a o ~ 10~ 2 ; we shall 
take an upper limit of go/o-o < 1/50. If the fluctuations 
in the Newtonian gravitational potential (which in linear 
theory are constant in time in an Einstein-de Sitter cos- 
mology) do not shrink, then go(z)/ao <, 10~ 2 (1 + z) 1 / 2 . 
This is an upper limit, and the accelerations could be 

2 This sortjpf behavior is not unique to the Lagrangian formulation; 
given by Eq. hi 
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significantly smaller at z = 3 in a MOND universe if 
the potential fluctuations have grown considerably since 
then. From Eq. || above, we see that an external field 
go J> ao/170 is required to modify the absorber internal dy- 
namics. Thus we can consistently consider external fields 
with 25 <, a /go <, 170. 

We have computed the total gas density in Lya ab- 
sorbers in both the isolated and non-isolated cases using 
the data of Hu et al. (1995) and Petitjean et al. (1993). 
In the isolated case we find f2 gas « 0.008 — 0.009. This is 
somewhat smaller than the il gas « 0.045 inferred from the 
same data in the CDM picture or using nucleosynthesis, 
but is still plausible - even if we assume that most baryons 
must be in gas in the IGM - considering the uncertainty 
inherent in the analysis and the neglect of undcrdense and 
collisionally ionized hot gas. In the case of a significant 
external field, we find 

n eas a 0.028 (^Jp) ~ 1/3 V. (12) 

giving 0.005 <J r2 gas <J 0.03; the higher end is in comfort- 
able agreement with observational constraints. 

For isolated absorbers, MOND also predicts that the 
low column density (~ 10 14 cm -2 ) absorbers have a rather 
small characteristic size: ~ lOkpc versus ~ 100 kpc for 
Newtonian gravity with CDM. We have argued that in a 
plausible external field, the absorbers should have sizes of 
20 — 80 kpc in their long direction, and half this in their 
short direction. Spectra of lensed quasars and close quasar 
pairs can be used to constrain the characteristic transverse 
sizes of Lya absorbers. On very small (few kpc) scales, ab- 
sorbers are virtually identical in both sight-lines (Smette 
et al. 1992; Dolan et al. 2000). On intermediate (tens 
of kpc) scales spectra are very similar but are currently 
too low-resolution to conclusively constrain the absorber 
properties across the sight-lines (e.g., Bechtold et al. 1994; 
Smette et al. 1995). On the largest scales, statistical anal- 
yses of the probability of detecting an absorber in both 
sight-lines lead to estimates of absorber 'sizes' of several 
hundred kpc (e.g., Dinshaw et al. 1995; Smette et al. 
1995; Crotts & Fang 1998; D'Odorico et al. 1998). These 
observations are slightly at odds with the MOND predic- 
tions unless go is near its upper limit, but some additional 
considerations must be kept in mind: 

1. Except on very small scales, the observations do not 
currently rigorously distinguish between correlated 
absorbers and a single absorber spanning two sight- 
lines. So we cannot directly compare the observa- 
tionally deduced sizes (several hundred kpc) to the 
Jeans length of absorbers of the observed column 
density (Eqs. | and [ll] ). 

2. Strictly speaking, the relations given in Eqs. and [ll] 
connect observed column density to radial dimension 
whereas the observations probe the transverse extent 
(or correlation). However, for an 'isolated' absorbing 
region these dimensions are unlikely to be very dif- 
ferent because absorbers should vary transversely on 

similar result can be derived from the simple formulation of MOND 
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a scale comparable to the local Jeans length, which 
is in turn comparable to the radial extent over which 
the absorption occurs. In the case of a strong exter- 
nal field, the absorbers should be elongated in the 
field direction by a factor of two. 

3. The properties of absorbers derived here assume that 
the absorbers are not far from local dynamical equi- 
librium, which would not be true for underdense 
and/or very large (3> 100 kpc) absorbers with dy- 
namical times exceeding the Hubble time. For iso- 
lated absorbers, it can be shown that regions of this 
size could not give rise to sufficient absorption un- 
less their dynamical time is shorter than the Hubble 
time. To see this, note that requiring idyn > £w_for a 
fixed L gives an upper limit on density via Eq. |2|, and 
this in turn gives an upper limit to column density 
of 



N H i < 4.4 x 10 11 cm~ 2 L 



3 

100 



tn 



10 



17. 



rp-0.76 

±4 

Tl2 ' 



where Lioo = i/100kpc. This means that isolated 
~ 100 kpc absorbers of column density ~ 10 14 cm~ 2 
cannot consistently be far from local dynamical equi- 
librium. Absorbers in an external potential go can 
be out-of-equilibrium if their column density obeys 
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Thus absorbers in a strong external field may be only 
marginally in equilibrium. 

In summary, we find that if the dynamics of Lya ab- 
sorbers in MOND are dominated by their self-gravity, then 
they are significantly smaller than observations indicate. 
If, however, the clouds are immersed in a (constant) accel- 
eration field of magnitude go, then their sizes can be sub- 
stantially larger, and the basic properties of Lya absorbers 
approach those predicted by standard gravity (with CDM) 
as go — ► (ri g as/^dm)ao ~ 0.16dQ) where f2 gas and flam are 
the density parameters in absorbing gas and dark matter 
in the CDM model. The importance of external fields in 
low- acceleration systems is both a methodological barrier 
and a saving grace of MOND. The accurate description of 
such systems requires a (necessarily ab-initio) calculation 
of the large-scale density field surrounding them, which is 
in turn impossible to perform rigorously without a cosmo- 
logical formulation of MOND which treats its conceptual 
problems (in particular the issue of which accelerations, 
and with respect to what, should be 'counted'). On the 
other hand, external fields can ensure that in the limit 
of extremely low accelerations the properties of isolated 
systems in MOND will not deviate wildly from their New- 
tonian counterparts. 

3. MOND AND CLUSTERS 

Like Lya systems, clusters of galaxies are gas-dominated 
(in the absence of dark matter), well observed (this time 
via X-ray measurements), should be in local hydrostatic 
equilibrium, and are in the MOND acceleration regime. 
They also have internal accelerations larger than expected 



ambient acceleration fields (and comparable to those near 
galaxies). Thus strong predictions regarding their struc- 
ture in MOND can be made using relatively simple argu- 
ments, as follows. 

For a gaseous system of mass M in hydrostatic equilib- 
rium, M is close to the Jeans mass Mj. For pure gas with 
MONDian gravity, 
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using Eq. ^. For comparison, the observed mass- 
temperature relation of clusters found by Mohr, Math- 
iesen, & Evrard (1999) is 



M ICM ~ 4 - 3 X 1Ql2 M 



IkeV/ 



(14) 



where Mj 5 ^ is the enclosed mass of the X-ray gas at the 
radius within which its density is approximately 500 times 
the cosmic mean. While only an order-of-magnitude esti- 
mate, Eq. [t3]would at first seem to be a remarkable success 
of MOND, since the CDM picture has difficulty explain- 
ing the observed mass-temperature relation in detail (see, 
e.g., Mohr et al. 1999; Finoguenov, Reiprich & Bohringer 
2001 and references therein). The fact that MOND can 
roughly account (to within a factor of two) for the mass 
discrepancy in clusters has been pointed out by Sanders 
(1999). But a hint of trouble is suggested by the fact that 
Eq. [lj applies for a particular radius, whereas for isother- 
mal clusters Eq. [l^ does not. A closer look reveals that 
serious problems arise when the temperature profile T(r) 
predicted by MOND for a given density profile p(r) is com- 
pared to observations. 

Consider a spherical system such as a cluster in hy- 
drostatic equilibrium. Then the temperature and density 
obey 

-— = -a, (15) 
p ar 

where a is the magnitude of the radial gravitational accel- 
eration, and P = (kT/pm p )p is the pressure. For a <C ao 
in MOND, this can be rewritten as 



dlogp dlogT prrip 
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a log r a log r kl 
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(16) 



where M(r) is the enclosed mass. This equation immedi- 
ately implies that if p(r) and T(r) are power laws, then 
T(r) oc y/M(r), as per the Jeans argument. More gener- 
ally, an increasing M(r) implies that T(r) increases as long 
as the sum of a p = dlog p/dlogr and ar = dlogT/dlogr 
changes more slowly than M(r) 1 / 2 . Thus isothermal den- 
sity profiles in MOND tend to have a core which contains 
most of the mass (and in which the logarithmic derivatives 
change relatively quickly) , then fall off more quickly than 
r~ 3 (see Milgom 1984). We will show that this is not the 
case for the X-ray emitting gas that dominates the mass 
in observed clusters, yet clusters are observed to be nearly 
isothermal outside of their central regions. 
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This constitutes a grave challenge to MOND which can 
be demonstrated using the specific cases of the Virgo, 
Abell 2199, and Coma clusters, for which gas density, stel- 
lar mass, and temperature profiles are available in the lit- 
erature. The analysis is shown in Figs, [j], || and [|. The 
first panel of Fig. |l| shows the cluster gas density depro- 
jected from a ROSAT X-ray emission profile by Nulsen & 
Bohringcr (1995), normalized to the critical density, as- 
suming a distance to M87 of 16Mpc. The density can 
be extended to r > 200 kpc using the fJ— model fit by 
Schindler et al. (1999) at large radii. The second (top- 
right) panel shows the integrated gas mass, as well as in- 
tegrated stellar mass from Schindler et al. (1999), with 
a contribution for M87 with M/Lb = 8 included from 
Giraud (1999). As discussed below, the details of these 
assumptions matter very little. The bottom-left panel 
shows the Newtonian acceleration at r, which is below the 
MOND parameter for r ^ 20 kpc; the MONDified accel- 
eration 3 is also plotted, along with the dynamical time 
y^l/a in units of the z — Hubble time. Given this infor- 
mation, Eq. [l6]can be used to predict T(r) given a starting 
T(tq). The lower-right panel shows this prediction, inte- 
grating inward 4 starting at tq w IMpc, with T(ro) taking 
values between l/10th and ten times the measured ASCA 
temperature there. If MOND were correct, one of these 
profiles should roughly match the observed temperatures, 
but none of them do. Figures || and || shows the same 
analysis for Abell 2199 and Coma, somewhat richer and 
more relaxed clusters. The results are quite similar. 

We have verified that the results are robust to reasonable 
changes in the distances to the clusters, the mean molec- 
ular weight, and the MOND interpolation formula used. 
We have also experimented with different profiles and nor- 
malizations for the stellar mass distribution; these do not 
significantly affect the results unless the stellar mass is so 
large [M/L 3> 20) as to imply the presence of dark mat- 
ter (dark matter is discussed below). Significantly larger 
values of the MOND constant (ao ^ 5 x 10~ 8 cms -1 ) help 
improve the fit because (as demonstrated by the agree- 
ment of Eqs. ^|| and [u]), the MOND 'Jeans Temperature' 
t mond _ ( Aimp /fc)(a GM) 1 / 2 roughly agrees with ob- 
served cluster temperatures at large radii. Increasing eto 
moves this agreement to intermediate radii, but still can- 
not yield a reasonable fit of the entire profile (and would 
be incompatible with the value required by galaxy rotation 
curves and make Lya absorbers even smaller). 

As alluded to above, the difficulty in accounting for the 
cluster data in MOND - even when ao is allowed to vary 
- can be understood in more general terms using an in- 
equality derived from Eq. [l6| which applies to any range 
[rijTa] over which the temperature is non-increasing: 



p. 165), giving 
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Convective stability requires that \o>t\ < fl^pl lest entropy 
gradients be erased on a sound crossing time (Sarazin 1988, 

3 We use the same interpolation formula in the transition region 
a[i(a/ao) = ajv, where = x/(l + x 2 ) 1 ' 2 . 
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In Virgo, for example, M(r) increases by a factor of « 25 
between 100 and 1000 kpc, where the gas is observed to 
be roughly isothermal. But then a p (r = 100 kpc) — 1.3 
requires (in MOND) that a p (r = 1000 kpc) w -3.9, while 
ol p ps —1.5 is observed. Both power-law indices and to- 
tal baryonic mass are very well constrained quantities, so 
this is a serious violation. Cluster gas density profiles are 
generally well-fit at large radii by /3-models of form 



P{r) = Pa 



3/3/2 



(19) 



In this case a(r)/a(ro) = 2/[l + (ro/r) 2 ] < 2. Therefore if 
M(r)/M(r ) > (10/3) 2 w 11 at any radius within which 
the cluster has a non-rising temperature, and within which 
the density profile is well fit by a j3— model, then MOND 
is violated. If T(r) is exactly constant, the constraint is 
stronger and M(r)/M(r ) > 4 violates MOND. 

Since the relevant properties of Virgo, A2199 and Coma 
(isothermal or radially declining temperatures, increas- 
ing M(r), and slowly changing d\ogp/d\ogr as in the (3- 
model) seem generic in clusters at large radii (Neumann 
& Arnaud 1999; Finoguenov, David, & Ponman 2000; Ir- 
win, Bregman, & Evrard 1999), it is very hard to see how 
to reconcile MOND with the observations. A few possi- 
bilities which do not seem able to satisfactorily effect this 
reconciliation are: 

Clusters are not in hydrostatic equilibrium: The z = 
Hubble time greatly exceeds the dynamical time (see 
bottom-left panel of each figure) and the sound-crossing 
time, inside ~ 1 Mpc, so hydrostatic equilibrium should 
hold within that radius (see Sarazin 1988 for some discus- 
sion). Moreover, simulations (albeit in the CDM picture) 
show that hydrostatic equilibrium and spherical symmetry 
are good assumptions in inferring the gravitational force 
given observed density and temperature profiles (Schindler 
1996). Finally, we note that the observed magnetic fields 
in clusters do not appear sufficiently strong to significantly 
affect the force balance or the inferred cluster bulk prop- 
erties (e.g., Goncalves & Friaca 1999; Dolag, Evrard, & 
Bartelmann 2001). 

Measured temperatures are incorrect: while tempera- 
ture determinations do have significant errors (especially 
in ROSAT data), the predicted MOND temperature pro- 
file disagrees by many a from the ASCA measurements 
in Virgo, from ASCA and BeppoSAX measurements in 
A2199, and from XMM measurements in Coma. More gen- 
erally, there is no indication in observations using ASCA 
(e.g., Markevitch et al. 1998; White 2000; Finoguenov, Ar- 
naud & David 2001), BeppoSAX (Irwin & Bregman 2000), 
or XMM (Arnaud et al. 2001a, b) that clusters have steeply 
rising temperatures at large radii. 

Efficient conduction causes clusters to be isothermal: 
Figure [l] shows gives the conduction timescale (Sarazin 

as Milgrom 1983b and Begeman, Broeils, & Sanders(1991), i.e. 



The integration can also be performed outward, matching the observed temperature at small radii, with essentially the same results. 
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Fig. 1. — Predicted MOND temperature profile for the Virgo cluster. Top left: X-ray gas density profile in units of the critical density, 
from Nulsen & Bohringcr (1995) (pluses). The dotted line is a Hernquist model fit at (and used here at) r < 200 kpc by Giraud (1999); the 

solid line is a /3-model, fit by Schindler et al. (1999) at large radii (solid line) and used here at r > 200 kpc. Top right: integrated mass in gas 
and galaxies. The mass of M87 is from Giraud (1999) with M/Lg = 8; we also add a component representing galaxies from Schindler et al. 
(1999). Bottom left: The Newtonian and MONDian acceleration (in units of 10 _9 cms~ 2 ) at each radius, showing that the cluster is deep 
in the MOND regime for r 2> 20 kpc. Also plotted are the MOND dynamical time tdyn an d the conduction timescale t con d, both in units of 
the 2 = Hubble time. Bottom right: Predicted temperature profiles, starting at 1 Mpc with temperatures between 0.1 and 10 times the 
(ASCA) observed temperature there. ROSAT (Nulsen & Bohringer 1995, deprojected) and ASCA (Shibata et al. 2001, projected) tempera- 
ture profiles with la error bars are shown for comparison. The dark, solid line is the MOND 'Jeans temperature', Tj = (/im p /fc)(aoGM) 1 / 2 . 
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Fig. 2. — Predicted MOND temperature profile for Abell 2199. Top left: X-ray gas density profile in units of the critical density, from 
Siddiqui, Stewart, & Johnstone (1998) (pluses; dotted line is a spline-fit), with the f}— profile of Markevitch et al. 1999 (solid line) fit at large 
radii and used here at r > 200 kpc. Adjusted to h = 0.7. Top right: integrated mass in gas and galaxies. The central galaxy NGC 6166 is 
taken (conservatively) to be a point mass with My = —23.47 (Gcbhardt et al. 1996) and M / Ly = 8; we also add a component representing 
galaxies proportional to the gas mass as per Cirimele et al. (1997) with amplitude chosen so to give total mass fraction in stars at Mpc 
equal to theirs. Bottom left: The Newtonian and MONDian acceleration at each radius, showing that the cluster is deep in the MOND 
regime for r ^> 30 kpc. Also plotted are the MOND dynamical time t^ yn and the conduction timescale t CO ndi both in units of the 2 = 
Hubble time. Bottom right: Predicted projected temperature profiles, starting at 2 Mpc with temperatures between 0.1 and 10 times the 
(ASCA) observed temperature there. BeppoSAX (Irwin & Brcgman 2000) and ASCA (Markevitch et al. 1999) temperature profiles with 2cr 
errors are shown for comparison. The dark, solid line is the MOND 'Jeans temperature', Tj = (/im p /fe)(aoGM) 1 / 2 . 
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Fig. 3. — Predicted MOND temperature profile for Coma. Top left: X-ray gas density profile (solid), using the 0— model fit of Mohr et al. 
(1999), adjusted to h = 0.7. The dashed line is the galaxy stellar mass density, using M/Ly = 8 and luminosity density from The & White 
(1988). Top right: integrated mass in gas and stars. Bottom left: The Newtonian and MONDian acceleration at each radius, showing that 
the cluster in the MOND regime for all radii. Also plotted are the MOND dynamical time t dyn and the conduction timcscalc t cont j, both in 
units of the 2 = Hubble time. Bottom right: Predicted projected temperature profiles, starting at 2 Mpc with temperatures between 0.1 
and 10 times the observed temperature there. The XMM temperature profile (Arnaud et al. 2001a) with 2a errors is shown for comparison. 
The dark, solid line is the MOND 'Jeans temperature', Tj = (/im p /fc)(aoGM) 1 / 2 . 
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1988) 



where 



5n e r k 
2k(T) ' 



k(T) = 2 x 10 



li 



IkeVy 



5/2 



ergs K 



(20) 



(21) 



This is a lower limit, because it neglects magnetic fields, 
which could increase the conduction time by a factor of 
between several and several thousand (see, e.g., Rosner & 
Tucker 1989 and Chandran & Cowley 1998, respectively). 
If icond <, ihub for some range in radius, we would expect to 
see a nearly isothermal temperature profile there. However 
(and regardless of whether or not £ conc i <; tdyn), the cluster 
density profile would respond to this conduction by read- 
justing to restore hydrostatic equilibrium on a dynamical 
time. But then (by an inverted version of the argument 
given above), MOND would predict a density profile with 
quickly-varying logarithmic derivatives, contradicting the 
observed density profile. In other words, Eq. 16 prescribes 



of dark matter. To do so, we have fit dark matter den- 
sity profiles of various (NFW, (3— model, Hernquist model) 
parametric forms 5 , while constraining the ratio ft of en- 
closed dark mass to enclosed gas mas at 1 Mpc. By low- 
ering this ratio we can find, for each parametric form, the 
lowest value of ft for which a plausible fit can be found. 
All three models yield fairly good fits if ft is left free, 
and yield ft « 4 - 5 for Coma, ft w 2 - 2.5 in A2199 
and ft w 1 — 2 in Virgo. When ft is fixed and progres- 
sively lowered, we find minimal allowed values of ft w 3 
for Coma, ft « 1.5 for A2199 and ft w 1 for Virgo. 6 The 
analogous ratios at 200 kpc take approximate minimal val- 
ues of 11, 13 and 7.5 respectively. The details of the fits 
depend upon h, oo, the MOND interpolation formula and 
(of course) the parametric form chosen, but the general re- 
sult that MOND requires dark matter of at least one and 
up to several times the gas mass within 1 Mpc, and about 
ten times the gas mass within 200 kpc, is robust unless 
h < 0.5 or a > 2 x 10 -8 cms~ 2 . 



a 1-to-l relation between p(r), M(r) and Tjr), demon- 
strated in Figs. |l| and [| Since the predicted MOND T(r) 
is not isothermal, an isothermal T(r) cannot match the 
observed p(r). 

Observed density profiles are incorrect: It has been pro- 
posed that a multiphase medium in clusters in which one 
component has a small filling factor could lead to an over- 
estimate of gas the density in X-ray measurements (e.g., 
White & Fabian 1995; Gunn & Thomas 1996). If this 
effect existed and were more severe at large radii, the ob- 
servationally inferred density profile could be flatter than 
the true profile, which could have more quickly-varying 
otp. However, studies of this effect find that the mass 
discrepancy is likely to be relatively small (<C 50%) and, 
moreover, less important at large radii (White & Fabian 
1995; White & Gunn & Thomas 1996; Nagai, Sulkanen, 
& Evrard 2000). The comparison between masses inferred 
from Sunyaev-Zel'dovich measurements (oc n e ) and X-ray 
emission measurements (oc n 2 ) also disfavor large correc- 
tions (Grego et al. 2001; Patel et al. 2000; Nagai et al. 

2000) . 

Based on observations of an apparent excess of extreme- 
UV emission, it has also been claimed that clusters may 
contain a large mass of warm (~ 10 6 K) gas (e.g., Lieu, 
Bonamente, & Mittaz 2000; Bonamente, Lieu, & Mittaz 

2001) . If this is the case, it would be extremely favor- 
able to the MOND hypothesis. Note, however, that an- 
other group, while finding EUV excess, finds a significantly 
smaller intensity (e.g., Bowyer, Korpela, & Berghofcr 
2001), and that FUSE observations place tight limits on 
warm gas in Coma and Virgo which are inconsistent with 
the proposed models (Dixon et al. 2001). 

There is dark matter in clusters: One might argue 
that there is a strongly concentrated component of bary- 
onic dark matter in clusters (Milgrom 1999; Sanders 
1999), leading to nearly- isothermal temperatures where 
that component dominates the gravitational mass. Us- 
ing our calculations we can estimate the required amount 

5 We have also tried several more general forms with no significant effect on the conclusions. 

6 Note that we are able to fit temperature profiles - in standard gravity - using NFW profiles with 7 < £i < 9 and concentration parameters 
between ~ 3 — 4 (for Coma) and ~ 7 — fO (A2199 and Virgo). The concentration for Coma is perhaps somewhat low but the results seem 
otherwise reasonable. 



The issues we have addressed in this section have been 
considered before by Gerbal et al. (1992, 1993) who com- 
puted M(r) using observed (but extrapolated) p{r) and 
the assumption of isothermality. This yielded M(r) some- 
what larger than the observed mass. Milgrom (1993) crit- 
icized this approach and suggested a different one similar 
to that employed here. Using the Coma cluster, The & 
White (1988) tested whether they could generate a set of 
p(r), M(r) and T(r) in MOND that could match the ob- 
servations available at that time. They succeeded, but 
only by using h = 0.5 and a — 2 x 10~ 8 cms~ 2 (which are 
not now observationally viable). We have reproduced their 
calculations, but find that even with their parameters, the 
now-available accurate temperature profile of Coma can- 
not be fit. 

Thus it seems that while MOND can account for the 
non-Keplerian form of galactic rotation curves (and - 
somewhat surprisingly - for the mass-temperature rela- 
tion in clusters), it cannot account for cluster density and 
temperature profiles in detail. The difference in success de- 
rives from the fact that - roughly speaking - in MOND an 
asymptotically isothermal temperature profile corresponds 
to a point mass, whereas clusters are extended yet still in 
the MOND acceleration regime. In CDM cosmology the 
isothermality of both galaxies and clusters is explained by 
assuming that both are embedded in an isothermal dark 
halo which dominates the mass. 

4. DISCUSSION AND CONCLUSIONS 

MOND as an alternative to the dark matter hypothe- 
sis has generally fared quite well when applied to galax- 
ies. However, galaxies can only test MOND in the limited 
regime of accelerations of (O.l-l)ao and physical scales of 
< 100 kpc. It is therefore important to discover whether 
MOND's success extends to systems at much lower ac- 
celerations and/or much greater physical scales. To this 
end, we have predicted various properties of Lya absorbers 
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and galaxy clusters in MOND, using hydrostatic equilib- 
rium arguments. Both types of systems are observationally 
well constrained, and well within the MOND acceleration 
regime. 

We find that as compared to their properties in the 
CDM picture, the Lya absorbers in MOND have some- 
what higher characteristic density and smaller character- 
istic size for a given column density, potentially in con- 
flict with absorption studies of quasar pairs and lenses. 
The magnitude of the effect, however, depends upon the 
(unknown) magnitude of the external acceleration field in 
which they embedded, since MOND violates the strong 
equivalence principle, 7 and if the ambient acceleration field 
is large then the predicted absorber properties can ap- 
proach those observed. When combined with the observed 
differential number density of absorption lines, the anal- 
ysis yields (for any assumed external field) a reasonable 
density of intergalactic gas (as does the CDM picture). 
A more accurate prediction of the detailed physical and 
statistical properties of Lya absorbers in MOND - which 
match observations in detail in the CDM case - will prob- 
ably require a (presently infeasible) ab-initio calculation of 
the large-scale acceleration field and the absorbers them- 
selves. 

We note also that the MOND Jeans mass depends only 
on temperature: 

Mj w 3.4 x 1O 6 M T|, (22) 

and gives a rather small Jeans mass in the 1GM which may 
dramatically alter galaxy formation in MOND (though we 
have not pursued that issue here). 

Simple arguments can also predict the mass- 
temperature relation in clusters, and the temperature 
profile of any single cluster given its observed gas density 
and enclosed mass profile. The MOND-predicted M-T 
relation is impressively close to the observed one, though 
this success seems coincidental, as it is sensitive to the 
radius at which the enclosed mass is measured (since in 
MOND the Jeans mass depends only on the temperature 
and the clusters are roughly isothermal). 



Stronger constraints can be derived using the observed 
mass, temperature, and density profiles of clusters: given 
hydrostatic equilibrium, MOND directly predicts the re- 
lation between the three quantities. For the form of mass 
and gas density profiles generally observed, MOND pre- 
dicts rising temperature profiles. In the specific cases of 
the Virgo, Abell 2199, and Coma clusters, we have shown 
that MOND's predictions strongly disagree with measure- 
ments from ASCA, ROSAT, BeppoSAX and XMM, and 
we see no reasonable way to effect a reconciliation without 
recourse to large amounts of (presumably baryonic) dark 
matter of an unknown type. For these clusters, the mass of 
such dark matter must exceed the gas mass within 1 Mpc 
by a factor of ~ 1 — 3, and by a factor of about ten within 
100 kpc. Moreover, we have argued that the discrepancy 
applies to clusters in general. This may be interpreted as 
a failure of MOND to describe cluster dynamics in terms 
of their observed baryonic content, or as a bold prediction 
(Sanders 1999) that we have so far observed only a mi- 
nority of their baryonic content. (Discovery of such dark 
matter would also constitute a serious crisis for the CDM 
model.) 

In conclusion, we find that although MOND can explain 
the rotation curves of galaxies in a simple and compelling 
way, it is less effective in extended systems such as clusters 
and (perhaps) intergalactic gas clouds, in which the visible 
mass cannot be described as a gravitational point mass 
when the system is in the MOND regime. This implies 
that dark matter (or perhaps some different modification 
of gravity) is required to accurately describe such systems. 
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